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Exploring an analytical expression for the convex roof of the pure state squared concurrence for 
rank 2 mixed states the entanglement of a system of three particles under decoherence is studied, 
using the monogamy inequality for mixed states and the residual entanglement obtained from it. 
The monogamy inequality is investigated both for the concurrence and the negativity in the case of 
local independent phase damping channel acting on generalized GHZ states of three particles and 
the local independent amplitude damping channel acting on generalized W state of three particles. 
It is shown that the bipartite entanglement between one qubit and the rest has a qualitative similar 
behavior to the entanglement between individual qubits, and that the residual entanglement in terms 
of the negativity cannot be a good entanglement measure for mixed states, since it can increase under 
local decoherence. 



I. INTRODUCTION 

The behavior of entanglement in systems of many par- 
ticles and the effects of decoherence on it is an important 
topic and far from being totally understood. If for one 
reason this is an issue in studying the transition from the 
microscopic to the classical world, it also has its practi- 
cal implications since the advent of Quantum Information 
and Computation; a scalable quantum computer is likely 
to need entanglement at a many particle level. 

A major obstacle in these studies is the concept of 
entanglement per si, as there is no established theory 
for multipartite entanglement (ME) and, even worse, no 
universal ME quantifier. One possibility is to analyse 
the entanglement between all possible bipartitions in the 
system (see [3], for example). But even in this case there 
is no easy path, given that most of bipartite entangle- 
ment measures do not have an analytical expression for 
mixed states of systems with even small Hilbert space 
dimension. The Entanglement of Formation and its re- 
lated quantity, the Concurrence, for example, are only 
exactly computable for 2x2 systems (two qubits) [2|. 
An exception is the negativity, which can be obtained in 
any dimension However for Hilbert space dimension 
higher than 2 x 3 it can underestimate the entanglement: 
it can be null when the state is entangled, defining what 
is called PPT entanglement or bound entanglement (two 
recent reviews on entanglement measures are [1, 0| ) ■ 

Somehow related to this approach to multipartite 
entanglement and relevant to it is the property of 
monogamy: if two qubits are maximally entangled, nei- 
ther of them can be entangled, in any way, with a third 
one (not even classically correlated 1 ). But it is in the case 
where the two particles are only partially entangled that 
the monogamy property may be more relevant to multi- 
partite entanglement. In particular, in the case of three 



qubits in a pure state Coffman, Kundu and Wootters @ 
obtained that C\, BC , > C\ B + C\ c , with Cab (Cac) 
being the concurrence of A with B(C), while C A (bc) ls 
the concurrence between A and BC, with the last treated 
as a single system. This inequality limits A's entangle- 
ment with B and C taken individually and do not allow 
both of them to increase indiscriminately, since all the 
quantities in the inequality varies from to 1. Besides es- 
tablishing this inequality Coffman, Kundu and Wootters 
also proposed to take the difference, C\^ BC ^—C\ B —C\ C , 
as a tripartite entanglement measure, which they called 
residual entanglement, t c abc . This monogamy inequality 
is also valid for the negativity and an analogous tripartite 
entanglement measure can be defined 

Even though for three particles the residual entangle- 
ment is well accepted as a tripartite entanglement mea- 
sure, most of the studies of the decoherence of multipar- 
tite entanglement only focused on the bip artite part of 
this entanglement @, @, M, ELI EI El (see [3] and 
[l6| for an exception and different approach). That hap- 
pens because in the case of mixed states despite the fact 
that the inequality and the residual entanglement can 
be established they involve a convex roof optimization of 
^a(bc) w ith no easy general solution. 

Here I make use of an analytical expression for the 
convex roof of the pure state C A ^ BC ^ for rank two mixed 
states, obtained by Osborne to study the monogamy 
inequality and the residual entanglement under decoher- 
ence for the first time. Both are investigated in terms of 
the concurrence and the negativity for generalized GHZ 
and W states of three qubits. I also compare the behavior 
of the N^ (BC) and its pure state convex roof extension. 



II. MONOGAMY AND MULTIPARTITE 
ENTANGLEMENT 



[1] By the other side a maximmaly classical correlation between A 
and B, also prohibits any of them of being entangled with a third 
particle. In PRA 69, 022309 the monogamy between classical and 
quantum correlation is investigated. 



One of the first investigations of the monogamy in- 
equality and its relation to multipartite entanglement was 
the work of Coffman, Kundu and Wootters in 2000 @|. 
There they first showed that for a pure state of three 



2 



qubits 



4 det p A >C 2 AB + C 2 AC , 



(1) 



with pa the reduced density matrix of the qubit A. To 
get the final inequality in terms of entanglement it is 
necessary to make use of two facts: 

1. the concurrence of a pure state of two qubits, A 
and B, is given by pa- 

2. even though the state space of AB is four dimension 
only two of these dimensions are needed to express 
a pure state of ABC (this comes from the Schmidt 
decomposition of the state ABC in the bipartition 
A|BC). 

Therefore A and BC can be treated as a pair of qubits 
in a pure state and 1\J det pa taken as the concurrence 
between A and BC, obtaining the monogamy inequality 



AB 



C 2 



AC- 



(2) 



This inequality allows one to define the residual entan- 
glement as 



'ABC — °A(BC) ~~ °AB 



r 2 

^AC- 



(3) 



As shown in their work, the residual entanglement does 
not depend on which qubit is chosen as A (the focus) 
and can thus be considered, as remarked by the authors, 
as a collective property of the three qubits, that measures 
an essential three-qubit entanglement. The GHZ state 
has t abc = C A ( BC j = 1, while a generalized W state, 
a|001) + /3|010) + 7|100), has t c A bc = 0. In this way it 
is usually said that the W state only contains bipartite 
entanglement, while the GHZ state only contains genuine 
tripartite entanglement. Note however that both states 
can be considered as genuine tripartite entangled as none 
of them can be written in an separable way, whatever 
bipartition is used. 

Before proceeding to the case of mixed states it is worth 
remarking some points which will be relevant. The con- 
currence of a pure state of two qubits is related to the 
linear entropy, Sl, of one of the sub-systems 2 



Cab = 2 Vdet p A = y/S L (pa) = ^2 (1 - tr [p 2 A ]). 

Another point that I should mention is that most of 
the difficulties in obtaining an entanglement measure for 
mixed states, comes from the fact that these are usually 
defined as convex roof extensions of a pure state mea- 
sure. Suppose, for example, that we have a well defined 



measure of entanglement for pure states, E(i(j). Then 
the entanglement of a mixed state p is defined as 

inf Vp,£W, 
{!>«,<#>»} ~ 

with the infimum taken over all possible decompositions 
of p in a mixture of pure states, p = ^ i Pi\4>i)(4'i\- The 
concurrence of a mixed state of two qubits, for exam- 
ple, is the convex roof of the pure state concurrence, or 
equivalently, the convex roof of the square root of the 
linear entropy, and in the especial case of two qubits an 
analytical expression can be obtained. 

That said, let me return to the monogamy inequality, 
but for mixed states. At first, when the three qubits 
are in mixed state Ca(bc) m n °t directly defined as all 
the four dimensions of BC might be used. Nonethe- 
less a convex roof extension of the squared concurrence, 

( C A(BC)) min = ^{p,^}T.iPrC 2 A{ BC) (&)> Can be USed 

to show that Q 
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Prom this inequality we can also define a residual entan- 
glement for mixed states of three qubits: 



T ABC ~ (Ca(BC)) 



C 2 



AB 



c 2 



AC- 



(5> 



However, in general, it will depend on which qubit was 
chosen as the focus. Thus for the case of mixed states 
it is necessary to use the average over the three possi- 
bles focus as the definition of residual entanglement, and 
that will be my procedure. From now on I will only use 



C 



A(BC) 



as the squared concurrence between A and BC, 
and when the state is mixed it is implicit that the convex 
roof extension was taken to obtain C A ^ BC ^ . 

The monogamy inequality can also be established for 
the negativity. As show by Ou and Fan in 2007 Q, for 
pure states of three qubits 



N 



A{BC) 



> NX B + N 



AC- 



(C) 



In the case of mixed states, despite N A ^ BC ^ being well 



definied without a convex roof procedure 3 , one has to 
mm to obtain, 



use ( N A(BC) 



(N 2 MBC) ) mm > N 2 AB + N 



'A(BC)/ ^ JV AB T ly AC- ( 7 ) 

Likewise the case of the concurrence Eq. [7] allows the 
definition of a residual entanglement in terms of the neg- 
ativity: 



'ABC 



(N 2 MB c)) 



^AB 



?AC- 
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[2] Some of the attempts to generalize the concept of concurrence to 
higher state space dimensions seem to preserve such equivalence 
I think one could even say that this equivalence can be used 
to such a generalization. Note that for pure states Sl(pa) is a 
valid measure of entanglement between A and B. 



[3] In Lee et al. [H] the convex roof of the negativity is explored as 
an entanglement measure. This is equal to the concurrence for 
two qubits, given the equivalence of the concurrence and nega- 
tivity for pure states with Schmidt rank 2, and can be seen as a 
generalized concurrence for higher dimensional systems. 
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Contrary to the case of the concurrence this residual en- 
tanglement depends on which qubit is chosen as A, even 
for pure states. Again, there is the possibility of tak- 
ing the average over the three possibles choices for the 
focus as the tripartite entanglement measure and, from 
now on, that is what I mean by r ABC . And it can be 
shown that this average is an entanglement monotone [3] 
for pure states. Since for pure states with Schmidt rank 
2, which is our case, the negativity is equivalent to the 
concurrence, then {N 2 A{BC) ) rnin = {C 2 A(BC) ) mm and I re- 
fer to { N \[BC)) mm as ^a(bc)- Contrary to the case of 
concurrence, the residual entanglement in terms of the 
negativity, t abc , is not null for the W state. Latter I 
will show that t abc is not a good entanglement quan- 
tifier for mixed states, since it can increase under the 
action of local operations. 

The difEcultty of studying the residual entanglement 
for mixed states, those that originate from decoherence, 
is to obtain an analytical expression for C\, BC y For- 
tunately, in 2002, Osborne obtained such analytical ex- 
pression for rank two mixed states [13]. Note that the 
convex roof of the squared concurrence, in general, is 
not equal to the square of the convex roof of the concur- 
rence: (C 2 ) min ^ ((C)" 1 ™) 2 . In fact as the pure state 
squared concurrence is a convex function then (tj 2 ) mm > 

((C) mm ) (the same is also true for the negativity). In 
the case of two qubits the equality was also shown by 
Osborne. Here I use the expression obtained by Osborne 
to study the residual entanglement in terms of the neg- 
ativity and concurrence of some states obtained after a 
decoherence channel has acted. I do not provide the de- 
tails of how to obtain this convex roof neither its final 
expression, referring the reader to the original article. 

It is important to remark that another possible defi- 
nition for the residual entanglement is to use the convex 
roof extension of the pure state residual tangle. This is 
different from the one employed here, has the advantage 
of being independent of which qubit is chosen as A, and 
was analyzed for a mixture of GHZ and W states with 
interesting results in [2p| ]. 

For the sake of completeness let me remark some other 
points about the monogamy inequality, before moving to 
the decoherence models I use. It is pleasant in the sense 
that it seems to reveal an entanglement structure: the 
entanglement of A with BC can be manifested in the 
form of bipartite entanglement with B and C taken in- 
dividually, plus an essential three-way entanglement en- 
volving all the three qubits. Besides this, it can be gen- 
eralized for a system of N qubits [2l|. Nonetheless, the 
inequality is not obeyed by all entanglement measures, 
being a negative example the Entanglement of Formation 
[f|, neither by higher dimensional systems [13] ■ In an y 
case, there are some attempts to explore this inequality 
and derivations from it to establish multipartite entan- 
glement measures for system of more than three particles 
M 0, EE H3, HI, Hi, with some progress and 
some drawbacks. Nowadays, in my opinion, it is not 



totally clear the real power and relevance of these pro- 
posals. 



III. DECOHERENCE MODELS 

It would be desirable to investigate the three paradig- 
matic types of decoherence: depolarization, phase damp- 
ing and amplitude damping. Of course, I also would 
like to apply these channels to general three (actually N) 
qubits states. But I have to restrict to the cases where 
the decohered state is rank 2, since this is the only case 
where an analytical expression for C A ^ BC ^ of mixed states 
is known. One possibility would be to use mixtures of 
generalized GHZ and W states, as they all have rank 2. 
The problem is that even these simple states, when suffer 
the influence of any of these three decoherence channels 
locally (each qubit is coupled to an independent channel) 
have their rank increased (note that a general mixture of 
three qubits can have rank 2 3 ). Even in the case of a 
pure GHZ or W state, the only exception is the local 
phase damping channel for the generalized GHZ state 
and the local amplitude damping for the generalized W 
state, which take their rank 1, as they are pure states, 
and changes it to 2 (see Tab. I). 





GHZ 


W 


DEP 


>2 


>2 


PD 


2 


: 2 


AD 


: 2 


2 



Table I: Decoherence Models and the rank of the mixed state 
generated when these are applied to generalized GHZ and 
generalized W pure states. DEP stands for Depolarization, 
PD for phase damping, and AD for amplitude damping. Here 
every particle is acted localy by a independent channel. Only 
the two cases where the rank is 2 are studied, as this is the 
only situation where an analytical expression for C 2 A ^ BC ^ is 
know 

I study the tripartite and bipartite entanglement in 
the two exceptional cases mentioned above 4 using the 
monogamy relation. Unfortunately, even though the 
GHZ and W states are paradigmatic tripartite entan- 
gled states, their distribution of entanglement is somehow 
trivial: the GHZ does not contain any bipartite entangle- 
ment between the individual parts, so that its tripartite 
entanglement is equal to the entanglement between A and 
BC. In contrast the W state only contains bipartite en- 
tanglement in the sense that the inequality is saturated 



[4] Another possibility is to apply the decoherence channel to only 
one of the qubits, but already in this much simpler situation 
just the amplitude damping channel when applied to generalized 
GHZ and W states generates rank 2 mixtures. However this last 
example will not be shown here, since it does not give qualitative 
difererences with the presented case. 
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and the residual entanglement is null (using the nega- 
tivity the inequality is not saturated for the W state). 
Because the decoherence channels only act locally and 
in an independent way, they are not able to create en- 
tanglement between the individual parts. Thus, the only 
possibility of having a significant different behaviour for 
the tripartite entanglement is to have rather different de- 
cays for the differents bipartitions. 

Decoherence and dissipation emerge when we couple 
the system of interest with another external system (usu- 
ally called environment, reservoir, bath or even ancilla) 
and only analyse the dynamics of the system of interest, 
ignoring the degrees or freedom of the external system. 
Mathematically these degrees of freedom are traced out 
and we end up with the reduced density matrix of the 
system of interest. The dynamics of this open system, 
which is influenced by the reservoir, can then be non- 
unitary and set in decoherence and/or dissipation. 

There are many formalisms to treat decoherence and 
dissipation in Quantum Mechanics. Here I will use the 
quantum operator formalism 5 , which describe the dy- 
namics experienced by the system as a completely pos- 
itive trace preserving (CPTP) linear map that acts on 
density operators: e(p). This map can include not only 
unitary evolutions, but also open ones and general mea- 
surements, besides other general operations. It turns 
out that any CPTP linear map can be written in an 
operator-sum representation: e (p) = ^ EipEj with 
J2i EiE} = I, being this last equality (a completeness 
relation) that guarantees the trace preserving property. 
The Ei are operators on the state space of the system of 
interest, usually called as Kraus operators, and can be 
obtained from the knowledge of the Hamiltonian govern- 
ing the whole system (interest + environment). 

The amplitude channel can represent a typical inter- 
action of a qubit with a zero-temperature reservoir in 
its ground state \0)r. In this interaction there is a finite 
probability, p, that the upper state |1) of the qubit decays 
to |0) and creates one excitation in the reservoir, which 
finishes at |1)_r. Of course there is a probability 1 — p 
that nothing happens. This interaction can represented 
by the map 



|l)|0) fl -> v / W|l)|0)ii + ^|0)|l) fi , 

and has the followings Kraus operators: E\ = 
(|0)(0| + vT=p|l)(l|) and E 2 = y/p\0)(l\. The prob- 
ability p is related to the interaction time and under 
Markovian approximation p = 1 — e~ r * with T a constant 
which depends on characteristics of the environment and 
its coupling with the system. 



The phase damping channel could also represent an 
interaction of the qubit with a reservoir. But, contrary 
to the amplitude damping, there is no loss of energy, but 
only information loss (in this sense the process is uniquely 
quantum). The eigenstates of the qubit do not change, 
but accumulate an unknown phase that destroys the rel- 
ative phase between them (loss of information). It rep- 
resents, for example, the elastic scattering between the 
qubit and the reservoir and is described by the following 
map 



10) 10)* 

11) 10)* 



^T^\i)\o) R + ^p\i)\i) R 



Note that the eigenstates of the qubit are not changed, 
but become entangled with the reservoir, and this will 
cause the loss of coherence in the qubit when the reservoir 
is traced out. The Kraus operators for the phase damping 
are: E 1 = vT=p(|0)(0| + |1)(1|), E 2 = ^(|0)(0|) and 
E3 = -y/p (|1)(1|). In the situation where the channel acts 
independently on each qubit of the system the formalism 
extends directly. 



IV. RESULTS 

First I present the behaviour of a generalized GHZ 
state under the action of three independent phase damp- 
ing channels on each qubit. In Fig. 1 the behavior of the 
entanglement of A with BC given in terms of the nega- 
tivity Na(bc) (dashed curve) and of the convex roof of 
the squared negativity C\^ BC ^ (solid curve) can be seen 
as a function of p: the convex roof decays faster than 
the negativity. However it may be more fair to com- 
pare C^( BC ) with the squared negativity N^ BC y This 
is not shown here, as the second becomes equal to the 
first (see footnote 3). Note that here the bipartite en- 
tanglement given by the convex roof is equal to the tri- 
partite entanglement given by the residual entanglement: 
t abc = t abc = ^A(BC)- F rom the graphics on the right 
of Fig.l it can also be checked that any generalized GHZ 
has less residual entanglement than the GHZ itself for 
any value of p and that the residual entanglement always 
decrease with p 

The second possibility is to look at the action of the 
amplitude damping in a generalized W state. In this case 



AT 2 



7^ ^A(BC)i but they have a qualitative similar 



[5] See chapter 8 of [3y| for more details and limitations. Recently 
results on the generality of completely positive maps as descrip- 
tions of dynamical process have been obtained in |'i ll |. 



behaviour. As an illustration, Fig. 2 shows both of them 
for the W state (upper curves in red) and a generalized 
one (lower curves in blue), which without normalization 
is |001) + 5|010) + 10|100). Note that the difference be- 
tween the two measures increases with decoherence at the 
beginning and then goes to zero again. In fact if N^ BC ^ 

is equal to zero, then C\^ BC ^ should also be, since there 
are no bound entangled states of rank 2 [32j . 

After viewing the difference between C^/ BC <) anc ^ 
^A(BC) I investigate the distribution of the bipartite en- 
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Figure 1: (color online) On the left N A (bc) (clashed) and 
^a(bc) (solid) for a state obtained when a phase damping 
channel acts locally and independently on each of the three 
qubits of an initial GHZ state (upper curves) and a generalized 
one, which without normalization is 0.2|000) + 1 1 1 1) . On the 
right C\( BC } for a family of generalized GHZ states, 1|000) + 
/3 1 1 1 1) . For the two states used on the left and the ones on 




P 



Figure 2: (color online) C 2 A( ^ BC ^ (solid) and JV|( BC ) (dashed) 
obtained when an amplitude damping channel acts locally 
and independently on each qubit of an initial W state (upper 
curves) and a generalized one, which without normalization 
is 001) +5|010) + 10|100). 

tanglement, using the concurrence as the measure of en- 
tanglement between qubits. In Fig. 3 I show this distri- 
bution in terms of the concurrence for the W state (left) 
and a generalized one (right), which without normaliza- 
tion is |001) + 2|010) + 3|100). There the behavior of the 
entanglement between one qubit and the rest for differ- 
ents focus, C\, BC y for example, under decoherence can 
be seen as solid curves. The entanglement between indi- 
vidual qubits is shown in dotted curves, while the residual 
entanglement is not show as it is null at the begginig and 
for any value of p. This indicates that the sum of the 
bipartite entanglement between the qubits decays at the 
same rate of the entanglement between one qubit and the 
rest. Fig. 4 is analogous but in terms of the negativity: 
we have N\ B instead of C\ B . In this case the residual 
entanglement t^ bc (dashed curve) may increase under 
the action of the local decoherence, indicating that this 
cannot be a good entanglement quantifier. 

V. CONCLUSIONS 

In sum, I first investigated the behavior of the squared 
negativity and its pure state convex roof extension under 
some models of decoherence for the generalized GHZ and 
generalized W states. While they are equal for the gen- 
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Figure 3: (color online) The distribution of entanglement in 
terms of the concurrence as a function of p, when an am- 
plitude damping channel acts locally and independently on 
each qubit of an initial W state (left) and a generalized one 
(right), which without normalization is |001) + 2|010)+3|100): 
Ca(bc)i Cfi(AC) an d Cc(ab) are ploted as solid lines from top 
to the bottom, respectively. C\ B , C\c an d C B q are given 
by dotted lines from top to the bottom, respectively. The 
residual entanglement r ABC , is null for any value of p. 
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Figure 4: (color online) The distribution of entanglement in 
terms of the negativity as a function of p, when an amplitude 
damping channel acts locally and independently on each qubit 
of an initial W state (left) and a generalized one (right), which 
without normalization is 001) + 2|010) + 3|100): C 2 A(BC) , 
Cb(ac) an d Cc(yis) are ploted as solid lines from top to the 
bottom, respectively. N AB) N ac and N BC are given by dot- 
ted lines from top to the bottom, respectively. The residual 
entanglement r ABC is show in dashed. Note that it may in- 
crease, showing that it cannot be an entanglement monotone. 

eralized GHZ state, that is not true for the generalized 
W state. However in this last case they exhibit similar 
qualitative behavior decaying with the decoherence even 
tough their difference can increase. 

I also studied how the distribution of entanglement 
behaves under these models of decoherence making use 
of the monogamy inequality for mixed states in terms 
of the concurrence and the negativity. For this aim I 
compared the behavior of the entanglement between one 
qubit and the other two treated as a single system with 
the one between individuals qubits. For the generalized 
W state I found out cases where the residual entangle- 
ment in terms of the negativity can increase under local 
decoherence even when all the bipartite entanglement is 
decaying, showing that it can not be a good entanglement 
measure. 

In relation to the behaviour of multipartite entangle- 
ment under decoherence these results show no qualitative 
difference between the bipartite and multipartite entan- 
glement for generalized GHZ and W states. However it 
would be interesting to be able to study these monogamy 
relations for states with richer entanglement structure 
than the GHZ and the W states. Would it be possible 



to happen that the bipartite entanglement between indi- 
viduals qubits goes to zero while the multipartite one is 
still finite, and maybe only goes to zero asymptotically 
(sudden death of only bipartite entanglement)? We ex- 
pect the rate of decay of both parts to be equal, since the 
tripartite entanglement should not increase under local 
operations. How about non-local environments? 

These studies could also increase our understandig of 
multipartite entanglement per si. For know, we showed 
that the residual entanglement in terms of the negativity 
as used here can not be considered as a good entangle- 
ment measure. But, as mentioned before, one ca also 
define the residual entanglement for mixed states as the 
convex roof of the pure state residual entanglement [2o| . 
And there are many other proposals of generalizations 
for higher Hilbert space dimensions and more than three 
particles ( (H, M, M, M, US H3, M, , for example). 
Therefore, it would be interesting to investigate the be- 
havior of these others measures under decoherence. 
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